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A biaxial rotational technique is proposed to improve the neutron irradiation uniformity for a large sample, 
and the theoretical method is established to predict and design the main parameters. The technique used a de- 
vice to rotate the target sample around two perpendicular axes simultaneously. Numerical calculations found 
that the lowest common multiple of the two angular speeds should be large enough to improve the uniformity, 
and the minimal experimental time should be no less than 600s. For a three-dimensional sample with a size 
of 20cm x 12cm x 14cm, the maximal non-uniform neutron irradiation factor of the sample is mainly deter- 
mined by the distance between the center of the sample and of the point neutron source. It was computed to be 
less than 10% when the distance was no less than 34cm. Experiments were carried out on the CFBR-II reactor 
and the experimental results were in good accordance with the theoretical analysis. As a result, the theoretical 
conclusions given above are reasonable and of reference value for the design of future irradiation experiments. 
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I. INTRODUCTION 


For the electronic materials and devices used in the space or 
nuclear radiation environment, neutron irradiation effects are 
widely studied to provide the foundation for radiation harden- 
ing design [1—4]. In practice, some of the electronic devices 
are relatively large, and it is often desirable to irradiate these 
large samples uniformly at high neutron fluence levels. 

A fast burst reactor is often adopted as the neutron 
source [2, 5, 6]. However, the neutron flux distribution is 
generally not flat by itself [7, 8], and the non-uniform fluence 
may become unacceptable for the study of the irradiation 
effects of a large sample [9, 10]. A more uniform neutron 
field can be achieved at a longer distance away from the 
source [7], but it may result in dramatic reduction in neutron 
flux in most cases [9, 10]. Ideal multiple radiation sources 
may improve the irradiation uniformity [10], but there is no 
such device available yet. A possible solution to this dilem- 
ma is to conduct irradiation experiments in the core of an ad- 
vanced fast-neutron reactor, as these reactors generally have 
large central cavities. The SPR III in the USA has a 16.5 cm 
diameter central cavity for the radiation testing of electron- 
ic parts and systems [2]. In Russia, the BR-1 reactor has a 
10cm x 18cm central cavity [11], and the BARS-5 reactor 
has two cores with diameters of 9cm and 11.5 cm, respective- 
ly. 

When the room in the core of a reactor is not enough for 
a large sample, as is the situation for the Chinese Fast Burst 
Reactor-II (CFBR-ID [12], an auxiliary device is required to 
improve the irradiation uniformity at a certain distance away 
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from the neutron source. A general rotator can be used to 
improve the irradiation uniformity, but is still not satisfactory 
for a large three-dimensional sample. 

This study aims to provide a theoretical foundation for the 
design of experiments with a novel biaxial rotational device, 
which will provide an enhanced neutron irradiation unifor- 
mity for the CFBR-II fast burst reactor. The biaxial rota- 
tional device rotates a sample around two perpendicular ax- 
es simultaneously. Numerical calculations and experimental 
verifications were carried out, and these results are of refer- 
ence value for the design of irradiation experiments. 


Il. THEORETICAL MODELS 


To provide a basis for the prediction and design of irra- 
diation experiments, the effect of the biaxial rotation on the 
irradiation non-uniformity of a three-dimensional sample was 
analyzed. The model of irradiation is shown in Fig. 1. 


Fig. 1. (Color online) Model of the irradiation experiment. 


The irradiated target is assumed to be a three-dimensional 
sample with a size of 20cm x 12cm x 14cm, and the loca- 
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Point (1,1,1) with w1= 0.1047 and w2= 1 
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Point (0,1,1) with w1= 0.1047 and w2 = 1 
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(f) 
Point (1,0,0) with w1= 0.1047 and w2= 1 
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Fig. 2. (Color online) Typical trajectory patterns of a point located on the sample. 


tion of the sample center is assumed to be the origin of co- 
ordinates. The biaxial rotational device rotates the sample 
around two axes simultaneously. The horizontal axis initially 
coincides with the X-axis and the angular speed around it is 
wı(rad/s). The vertical axis always coincides with the z axis, 


x£ = zo cos(we X t) — yo cos(wy x t) 


y = xo sin(we x t) — yo cos(wı X t) c 


z = yo sin(wı x t) + zo cos(wı x t) 


Obviously, it would be difficult to treat a three-dimensional 
neutron source as the CFBR-II reactor. A point neutron 
source located at R = (0, — R, 0) is considered in this model, 
and both scattering and absorption of the neutron field due to 
the sample are also neglected. The neutron fluence rate for a 
point located at 7 = (x,y,z) is (T), which is proportional 
to =a As a result, the average neutron fluence rate over 
time T for a point initially located at (xo, yo, 20) is 


iy? ig? 4 
(Xo, Yo, Zo) = 7 odt x zf FR 


The average neutron fluence rate at the center of the sample 
is selected as the reference value, (0), and the non-uniform 
neutron irradiation factor of this point is defined as follows, 


_ [8(r) - 20) 
(0) 
The maximum ņ for any point on the irradiated sample is 


regarded as the 7) of the sample. The goal of theoretical anal- 
ysis is to get a relatively small 7, e.g. less than 10% of the 


(2) 


x 100%. (3) 


and the angular speed is wa(rad/s). 

Before analyzing the irradiation non-uniformity, the trajec- 
tory of a point on the sample must be computed. For a point 
initially located at (xo, yo, 20) on the sample, we computed 
the time-dependent space coordinates of the point as follows. 


sin(w2 x t) + zo sin(wı x t) sin(we x t) 


os(w2 X t) — zo sin(w1 x t) cos(we x t), 


(1) 


sample considered, and the distance away from the neutron 
source, R, should be as small as possible to provide a rela- 
tively high neutron fluence level. 

The theoretical expressions were solved numerically by 
Matlab software. 


I. NUMERICAL CALCULATIONS 


A. Trajectory patterns 


As the sample center coincides with the origin of coordi- 
nates, the trajectory of a point on the sample is always locat- 
ed on a spherical surface. If the trajectory is more uniform, 
the non-uniform factor 7 should be smaller. The trajectory 
depends on both initial locations and the ratio of w2/w (if 
w 1 # 0). Typical trajectory patterns are shown in Fig. 2. 

When the angular speed, w2, around the vertical axis is set 
to be 1 rad/s, the change in w, would result in different kinds 
of trajectory patterns as shown in Fig. 2(a)— 2(d). For most 
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cases, if the initial coordinates of a point on the sample are 
located neither on the Z axis nor on the middle horizontal 
plane, the trajectory pattern is similar to that of the point ini- 
tially located at (1, 1, 1). The trajectory patterns in Fig. 2 
and other unlisted results suggest that the distribution of the 
trajectory path becomes more uniform as the lowest common 
multiple of w; and w2 becomes larger. The trajectory line is 
sparsely distributed if the w; is equal to 1 rad/s, 0.2 rad/s, or 
0.1 rad/s, but a uniformly distributed barrel structure can be 
obtained if wı = 0.1047 rad/s. A more uniform trajectory 
pattern is supposed to result in a smaller 7, and this assump- 
tion will be confirmed in the following text. 

The trajectory patterns also depend on initial location. On 
the condition of w; = 0.1047 rad/s and wə = Irad/s, the 
point initially located on the yz plane would result in a spher- 
ical trajectory pattern with an increase in time, and the point 
initially located on the x axis would result in a simple two 
dimensional circle. 


B. Minimal experimental time 


Point (10,6,7) 


400 600 
Time (s) 


Point (10,0,0) | 


400 600 800 
Time (s) 


% 200 1000 


Fig. 3. (Color online) The 77 as a function of time. 


Based on the relationship between the angular speeds and 
the trajectory patterns, the angular speeds of the biaxial rota- 
tional equipment were set to be w; = 0.1047, w2 = 1. If 
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the distance away from the neutron source is 36 cm, the non- 
uniform neutron irradiation factor, 7, as a function of time for 
the points initially located at (10,6,7) (Coordinates are in cen- 
timeters, so as the following text) and (10,0,0) are shown in 
Figs. 3(a) and 3(b), respectively. It can be observed that the 
oscillation of 7 becomes moderate with an increase in time 
for both cases, and the maximum difference of 77 is less than 
0.4% after about 600 s. As a result, the minimal experimental 
time should be no less than 600 s under the condition given 
above, and this requirement is generally easy to meet. 


C. Minimal R 


—— Point (10,6,7) 
-Point (10,0,0) 


Fig. 4. (Color online) The 77 as a function of the distance away from 
the neutron source. 


R has the greatest influence on the neutron fluence levels 
for most of the irradiation experiments, and it should be as 
small as possible to save experimental time. Under the con- 
dition of wı = 0.1047 and wə = 1, the theoretic value of 
7 was shown in Fig. 4 as a function of R for the points ini- 
tially located at (10, 6,7) and (10, 0,0). In order to alleviate 
the deviation caused by time, the 7 provided is the average 
result around time 600s. As expected, the 7) decreases quick- 
ly with the increase of the distance away from the neutron 
source R for both cases. The 7 for point (10, 6,7) is less 
than that for point (10,0,0), and R should be no less than 
34cm if 7 < 10%. Since the CFBR-II is treated as a point 
neutron source, we conjectured that R = 36cm might satisfy 
n < 10% under experimental situation. 


D. Angular speeds of the biaxial rotational equipment 


When R = 36cm and w2 = 1, the 7 of point (10, 6,7) and 
point (10, 0,0) as a function of w, were computed, as shown 
in Fig. 5. It can be observed that w; has little effect on the 7 
of point (10, 0,0) because the corresponding trajectory path 
is always a circle, but the 7 of point (10, 6,7) is sensitive to 
the selection of w; due to the change of trajectory patterns as 
shown in Fig. 2. The results suggest that the distribution of the 
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Fig. 5. (Color online) The 77 as a function of wi when w2 = 1. 


trajectory path should be as uniform as possible to enhance 
neutron irradiation uniformity. In the design of the irradiation 
experiment, the lowest common multiple of wı and wz should 
be large enough, e.g., when wı = 0.1047 and w2 = 1. 


E. Spatial distribution of 7 


The above discussions generally focused on the 77 of typical 
points when wı = 0.1047, w2 = 1, and R = 36cm. Based 
on analysis of the trajectory pattern of a point initially locat- 
ed at (20, yo, zo), the spatial distribution of 7 was shown in 
Fig. 6(a) to evaluate the 7) of the sample under the same con- 
ditions, where ry, = y yf + zê. It can be observed that ryz 
has little effect on 7, suggesting that the 7) is almost the same 
for a sample with a larger yo and zọ. The maximum 7) mainly 
depends on the value of zp. More explicitly, the computed 7 
is less than 1% when || is less than 3.5 cm, and the 77 is less 
than 8.5% when |xo| is not larger than 10 cm. 

For the generally adopted single axis rotator (around axis 
Z), the typical space distributions of 7 were shown in Fig. 6(b) 
with wı = 0 and w = 1. It can be observed that 7) depends 
on both zo and rey (ryz = y yê + 2), and the maximum 7 
is about 12%. For the static state irradiation experiments, the 
spatial distribution of 7 on the horizontal middle plane was 
shown in Fig. 6(c). 7 is very sensitive on the space location, 
and has a maximum value of 44%. As a result, we demon- 
strated that a sample rotating around two perpendicular axes 
simultaneously can result in significant enhancement of the 
irradiation uniformity. 


IV. EXPERIMENTAL VERIFICATION 


The schematic diagram and a picture of the biaxial rota- 
tional equipment are shown in Fig. 7. The real-time opera- 
tional state of the equipment can be monitored and controlled 
by a computer, and the angular speeds can be easily changed 
according to experimental requirement. 

The experimental layout is designed to be in accordance 
with the theoretical model given in Fig. 1, and the distance 
between the centers of the sample and the CFBR-II reactor 
is 36cm. The neutron fluence was measured by the neutron 
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Fig. 6. (Color online) Spatial distributions of 7 when (a) wı = 
0.1047 and w2 = 1; (b) wi = 0 and w2 = 1; and (c) wi = w2 = 0. 


activation method and nickel discs were adopted. 28 nickel 
disc were stuck on the 7 faces (6 outer faces and 1 vertical 
inner face) as shown in Fig. 1. With wı = 0.1047 and w2 = 
1, the maximal 7 measured in the experiment was 9.3%, in 
accordance with the theoretical value of 8.5%. For the static 
state irradiation experiments with w; = w2 = 0, the maximal 
n measured in the experiment was 49.0%, which also agrees 
with the theoretical value of 44%. As a result, the theoretical 
analysis given above is reasonable and of reference value for 
the design of future irradiation experiments. 


V. CONCLUSION 


Because the room in the core of a reactor may be not e- 
nough for a large sample, an auxiliary equipment is required 
to improve the irradiation uniformity at a certain distance 
away from the neutron source. A biaxial rotational technique 
is proposed to improve the neutron irradiation uniformity for 
a large sample, and the theoretical method is established to 
predict and design the main parameters. The lowest common 
multiple of the two angular speeds should be large enough to 
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Fig. 7. (Color online) (a) The schematic representation and (b) the picture of the biaxial rotational equipment. 


improve the uniformity, and the minimal experimental time 
should be no less than 600s. For a three-dimensional sam- 
ple with a size of 20cm x 12cm x 14cm, the maximal non- 
uniform neutron irradiation factor of the sample is mainly de- 
termined by the the distance between the centers of the sam- 
ple and the point neutron source. It was computed to be less 
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